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ABSTRACT 



Robotic telescopes present the opportunity for the sparse temporal placement of observations when period searching. We address 
the best way to place a limited number of observations to cover the dynamic range of frequencies required by an observer. We 
show that an observation distribution geometrically spaced in time can minimise aliasing effects arising from sparse sampling, 
substantially improving signal detection quality. The base of the geometric series is however a critical factor in the overall 
success of this strategy. Further, we show that for such an optimal distribution observations may be reordered, as long as the 
distribution of spacings is preserved, with almost no loss of quality. This implies that optimal observing strategies can retain 
significant flexibility in the face of scheduling constraints, by providing scope for on-the-fly adaptation. Finally, we present 
optimal geometric samplings for a wide range of common observing scenarios, with an emphasis on practical application by the 
observer at the telescope. Such a sampling represents the best practical empirical solution to the undersampling problem that 
we are aware of. The technique has applications to robotic telescope and satellite observing strategies, where target acquisition 
overheads mean that a greater total target exposure time (and hence signal-to-noise) can often in practice be achieved by 
limiting the number of observations. 
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1. Introduction 

Recent years have seen the rise of research-gr ade robotic 
telesc opes such as the Liverpool Telescope l|Bode et all 
Il999|) . as well as a sharp increase in the total num- 
ber of autonomous telescopes active across the world 
(|Hessman et al.ll2004|) . The emergence of these telescopes 
and the flexibility of their operation jusitifies renewed in- 
terest in the optimum scheduling of observations. Our aim 
is to discuss in detail one such strategy and to demonstrate 
its optimal properties for our project. Individual robotic 
facilities offer new scheduling possibilities to the observer 
by relaxing the requirement that a human be physically 
present at the telescope at all times. Robotic telescopes 
offer the potential for ideal sampling of periodic signals, 
since observations may be requested more sparsely and 
over a longer timescale than a traditional continuous ob- 
serving run. 

Furthermore, when unmanned telescopes are linked by 
a communication infrastruct ure, new ways of perform- 
ing science become feasible ll Allan et alJ 20041) . For ex- 
ample, the diurnal sampling problems that plague single- 
site observations of periodic variables can be avoided if 
an observer can acquire data from a second, geograph- 



ically distant site. Continuous monitoring projects such 
as the Whole Earth Telescope l|Nather et alJll99fJ). and 
the B iSON network for solar observations l|Chaplin et all 
1996) indicate the success of this approach. 

Thus we arrive at a central question in the study of 
periodic stars: given access to round-the-clock observing, 
how can we best sample a potentially periodic star, whose 
luminosity oscillates with some unknown frequency, such 
that we can minimise the effects of aliasing, while pro- 
viding a suitably long baseline with which to acquire the 
periods of interest, all while maintaining the quality of our 
data? In particular, how should we choose when to observe 
in those situations where the total number of observations 
that may be made is a limited resource? This situation is 
commonly the case for robotic telescopes, where target 
acquisition overheads, which include the slew time and 
the CCD readout time, as well as scheduler requirements, 
can preclude large numbers of observations. Although this 
question reflects our interests in using robotic networks to 
study variable stars, it is equally applicable to satellite- 
based observations. 

The astrophysics of the problem limits the frequencies 
of interest to an astronomer, allowing limits to be placed 
on the range of frequencies which must be correctly re- 



2 



Eric S. Saunders et al.: Optimal placement of a limited number of observations for period searches 



constructed. Then for an evenly sampled dataset, the in- 
put signal may be correctly reconstructed if the sampling 
frequency is at least twice the highest frequency f max to 
be detected. Assuming the astronomer wishes to see two 
cycles of any modulation then the lowest frequency de- 
tectable is given by 2/T, where T is the duration of the 
total sampling interval (i.e. the run- length). The value 
of the required sampling frequency z/n for equally spaced 
data can in some way be viewed as the Nyquist frequency 
of the dataset, given by 

^max < N/2T = Z/n, (1) 

where N is the number of observations. The problem be- 
comes more interesting, however, for the case in which 
there are not sufficient datapoints to fulfill this condition. 
Then the question is whether, under these conditions, we 
can sample in such a way that we can recover a signal of 
acceptable quality. Alternatively, we can reverse the prob- 
lem: how many datapoints do we need in order to recover 
an acceptable signal? In what follows, we show that there 
is scope for significant optimisation of observation place- 
ment in this undersampled regime. 

2. Previous work 

In the field of time domain astronomy, there is a 
relative dearth of literature regarding the best way 
to sample data in time. In contrast, much attention 
has been focussed on the prob lem of signal extraction 
from fixed, finite datasets ( e.g. iHorne fc Baliunasl 1*1986: 
ISchwarzenberg-Czernvll2f)03h . This is perhaps unsurpris- 
ing. The vast majority of such datasets are acquired in 
the classical situation of a single extended observing run, 
where an astronomer physically travels to a remote ob- 
serving site, accumulating data continuously each night 
for the duration of the run. In this case the observer has 
relatively little choice about how to space the observa- 
tions, and normally opts for the safe strategy of using all 
of the available time for continuous observing, acquiring 
a l arge total numbe r of observations. 

iDeemind lll975|) showed that the observed Fourier 
transform of a discrete dataset can be expressed as the 
convolution of the true Fourier transform F^(u) with a 
spectral window 6^{y) that fully describes the interfer- 
ence effects between the sampled frequencies, such that 

Fs(y) = F(v) * (2) 

where F(u) is the complex Fourier transform of a function 
fit) and is defined as 

/+00 
f(t)e i2 ^dt, (3) 
-00 

and the spectral window is given by 

JV 

&{y) = Y,^" n \ (4) 
fe=i 



where N is the total number of observations and tk is 
the set of observation times. The window function is even, 
so only positive frequencies are evaluated. The physical 
interpretation of the window function can be seen by con- 
sidering the Fourier transform of an infinite sine function 
(a delta function) convolved with the window function, 
leading to the observed Fourier transform. Adjusting the 
set of observation times tk al ters the window function, and 
hence the Fourier transform. IDeemind l|l975|) showed that 
aliasing effects in the observed Fourier transform arising 
from the window function could be mitigated by intro- 
ducing irregularities into the sampling pattern. By careful 
choice of sampling it is possible to reconstruct higher fre- 
quencies than i^n, the highest frequency that could be ac- 
curately reconstructed if the same number of observations 
were equally spaced. 

IDeemind l|l975|) illustrated the behaviour of the win- 
dow function through the use of a simple single-variable 
model for the N = 25 case, empirically deriving a form 
for the spacing based on v isual inspection of t he resulting 
spectral window function. iMould et al.l l|200(t) adopted a 
pseudo-geometric spacing as a way to maximise the uni- 
formity of their phase spacing, but no details of the im- 
plemented optimisation strategy were provided. 

The Rayleigh resolution criterion may be written as 
AvAT = 1. It means that observations of length AT 
yield the frequency of an oscillation with an error of 
Av — 1/AT, at worst. After an additional time interval 
oAT this yields an error in the cycle count AvaAT = a. 
Thus to keep the cycle count error less than a, subse- 
quent observations ought to be scheduled no later than 
(l + a)AT from the start. Recurrent and rigid application 
of this known principle yields a geometric series schedul- 
ing (we are grateful to the referee for suggesting this line 
of argument). 

We shall show that the choice of geometric base is cru- 
cial for such a strategy to succeed. The general question of 
what form an optimal uneven spacing should take, when 
the number of observations and the length of the observ- 
ing run are free parameters remains, to our knowledge, 
unanswered. 

Here we present a detailed analysis of the application 
of a geometric spacing technique to the undersampling 
problem. Although our sampling is not an analytic so- 
lution to the general problem of arbitrary pla cement of 
the da taset tk, it outperforms both the model of lDeeminel 
and pure random distributions, and thus represents 
the best current heuristic solution to the sampling prob- 
lem. In Section |3 we present the sampling strategy and 
quality metrics by which it is judged. In Section Elwe com- 
pare t he geometric sampling to the sampling of IDeemind 
(197,3) an d to random sampling and show that noise has 
no impact on the sampling strategy. We also show that 
the window function retains similar characteristics under 
random reordering of the spacings, an important result for 
practical scheduling scenarios. In Section we present a 
graph indicating optimal geometric samplings for observ- 
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ing runs spanning the range from 10 to 500 observations. 
We present our final conclusions in Section 



3. Method 

We generate our observation times with a geometric dis- 
tribution using the scheme 



1 



T. 



(5) 



where tk is the (time) value of the point in the series, x is 
the base of the series, k is the number of the current data 
point, N is the total number of observations and T is the 
duration of the total observing interval. This produces a 
geometrically spaced distribution such that < tk < T. 

The parameter x, which we call the geometric base, 
may be arbitrarily chosen. Since the wi ndow function is 
in general complex, for comparison with iDeemin 
we follow the convention of examining the normalised am- 
plitude A(v), defined as the square of the amplitude of the 
window function, normalised by the square of the number 
of datapoints, such that 



A{y) 



N 2 



(6) 



We define the optimal base, x opt , as the value for which the 
overall degree of structure present in A{v) is minimised. 
We use the root mean square (RMS) of A{v) as a straight- 
forward metric for identifying structure. Minimising the 
RMS thus represents the desire of the astronomer to 
achieve relatively even sensitivity to periods. 

In order to retain sensitivity to the higher frequency 
structure we evaluate A{y) for frequencies in the range 
0.1 < i^sr < 5 (in units of the Nyquist frequency for an 
equivalent even sampling). In practice, the range of the 
window function that needs to be considered will vary de- 
pending on the degree of undersampling of a given dataset 
(which depends on the number of observations and the 
length of the dataset), and the maximum frequency to be 
detected. Since what matters is the frequency with respect 
to the total time interval, we define the dimensionless rel- 
ative frequency as 



vT. 



(7) 



Then the relative limiting frequency in units of the relative 
Nyquist frequency is simply v max /v^. 

4. Results 

4.1. The best sampling 

generated 25 datapoints according to the 

formula 



fc- 1 /" (fc=1...12) 
(25 -fc)- 1 /" (fc = 13 . . . 24). 
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Frequency (Nyquist units) = v/v v 



Fig. 1. The window function (where the amplitude is 
given by \8(y)\ 2 /N 2 ) produced by Deeming's distribution, 
with a close to optimal value of a = 1.0. The RMS is 
0.0361. 

Frequency times Run Length = i/ rel 
20 30 40 




(8) 
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Fig. 2. The window function for a typical random distri- 
bution of 25 data points. The RMS is 0.0386. 



Figure ^ is the window function for Deeming's distribu- 
tion with a = 1.0 - a near-optimal value by visual inspec- 
tion and RMS comparison. Figure [21 presents the window 
function for a randomly spaced set of 25 datapoints. The 
typical size of the interference structures at high frequency 
are comparable, as indicated by the similar RMS values 
of the two graphs. Figure 01 is the window function for 
an optimal geometric spacing, where the value of x was 
chosen according to the minimum of Figure 01 The geo- 
metric spacing exhibits much better high frequency sup- 
pression, with an RMS roughly half that of either of the 
other spacings, at the price of a slightly broader true peak 
(centred by construction at v = 0). This broadening is a 
consequence of the undersampling. Although we can push 
the effective limit for frequency reconstruction to higher 
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Fig. 3. The window function for the optimum geometric 
sampling of 25 data points, for x = 1.124. The RMS is 
0.0201. 
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Fig. 4. RMS as a function of geometric base, for N = 25. 
The optimum value is x = 1.124. 



frequencies, we sacrifice some of our knowledge of lower 
frequencies in order to do so. In particular, the precise 
frequency of the true signal is less clearly represented in 
the window function. 

Note that the choice of geometric base is critical to the 
success of the sampling. Figure [5] shows how structure can 
be reintroduced into the window function by a poor choice 
of geometric base. This situation occurs because a sub- 
optimal geometric base essentially 'wastes' information, 
oversampling some frequencies and undersampling others. 
In the pathological case of even sampling, every sampled 
frequency is a multiple of the true frequency, leading to 
massive redundancy. For datasets with more than 25 dat- 
apoints, the correct choice of geometric base is much more 
important. Figure HO plots RMS as a function of geomet- 
ric base for datasets with 25, 100 and 500 observations. 
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Fig. 5. Window functions for three geometric spacings of 
25 datapoints. The top panel illustrates aliasing effects at 
a near linear sampling of x = 1.0094, while at x = 1.3 
(bottom panel) the effects of near frequency interference 
are apparent. The optimum geometric spacing for this 
number of datapoints is shown for comparison (x = 1.124, 
centre panel). 
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Fig. 6. Comparison of the variation of RMS with geomet- 
ric base, for three dataset sizes. When n — 25, a relatively 
broad range of bases lie close to the minimum RMS. For 
larger numbers of datapoints, the choice of optimal base 
is much more tightly constrained. 



The larger datasets possess much narrower minima 'val- 
leys' than the n — 25 dataset. At these dataset sizes, any 
substantial deviation from the optimal base will lead to 
rapid degradation of the window function. 
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4.2. Noise 

Whilst in the idealised case without noise a suitable geo- 
metric sampling is the correct sampling, we must consider 
the effects of noise on the result. If noise dominates the 
periodogram, then there may be little point in optimis- 
ing in the manner described. A range of equally plausible 
'optimal' geometric bases will exist, since the quality of 
the final observed Fourier transform is noise- limited rather 
than window-function limited. 

When making actual observations, an observer must 
decide how to divide a fixed total exposure time. This 
choice impacts signal-to-noise, since for most practical 
observing scenarios the signal-to-noise increases as the 
square root of the total exposure time. However, the win- 
dow function has no noise dependency. The quality of a 
specific window function as a function of frequency is de- 
termined entirely by the choice of sampling time. This 
means that noise will only be of practical concern to the 
observing str ategy if i t is gr eater than the mean ampli- 
tude of what DccmingJ l|197q) refers to as 'near-frequency 
interference', and which we shall term spectral leakage. 

To demonstrate the practical effects of noise on our 
ability to recover a signal from a set of observations, we 
generate a set of artificial lightcurves, each a sinusoid of a 
single frequency, that span the range from I- 1 00 cycles per 
observing run. We define the relative period as the ratio 
of the period P to the observing run T, such that 



.Pre 



P 

T 



1 

frel' 



(9) 



Thus, we consider the range 0.01 < P rc i < 1. For 
each gene rated period, we calculate the Scargle peri- 
odogram llLombl Il97rlt IScardei Il982ft . as formulated by 
iHorne k Baliunasl l|l986|) . 

The quality of the periodogram is assessed by calcu- 
lating the peak ratio, A, which we define to be the ratio of 
the powers of the peak corresponding to the true period 
and the highest other peak in the periodogram. There are 
three regimes of interest for this metric. Figure illus- 
trates the 'alias' regime, where the metric indicates the 
extent to which power has been lost to other strong sig- 
nal candidates. When aliases have been largely suppressed 
or when the overall level of noise in the periodogram is 
large, the metric instead describes the prominence of the 
peak relative to the spectral leakage, or to the background 
noise. These two situations are illustrated in Fig. IHlandEl 
respectively. The relationship to the alias property or peak 
amplitude leads us to denote this metric A. 

The accuracy to which a period detection is required 
changes A. An observer would not be interested in all 
nearby peaks or sub-peaks, but the scale on which struc- 
ture in the periodogram is important is obviously highly 
dependent on the nature of the science programme. For 
all our simulations, we consider periods within 10 per 
cent of the true period to be adequate, in line with typ- 
ical accuracy require ments for rotation rate histograms 
( Littlefair etai1l2005j) . We assume that our periods have 
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Fig. 7. Periodogram illustrating how the A metric can be 
applied as a measure of relative alias strength. A\ and 
A 2 are the amplitudes of the true peak (in this case, the 
highest peak) and the highest other peak in the transform. 
The ratio 4 1 is our A metric. 
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Fig. 8. Periodogram illustrating how the A metric, defined 
as the ratio can represent the power at the true period 
relative to the spectral leakage. No noise has been added 
to this signal; interference effects are due to the (relatively 
poor) choice of sampling. The signal is a simple sinusoid 
with P re i = 0.10163, sampled with 100 datapoints using a 
geometric base of x = 1.10, randomly respaced. The inset 
gives the amplitude of the corresponding window function, 
as a function of frequency 



an infinite coherence length, i.e. that the phase informa- 
tion does not change over time. 

For a 100 observation dataset, we apply noise equal 
to one third of the full amplitude of the signal (Figure 
HUfl and noise equal to the amplitude of the signal itself 
(Figure I12f) . This yields a signal-to-noise of 3.33 and 1 
respectively. It should be noted that this is not the signal- 
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Fig. 9. Periodogram illustrating how the A metric, de- 
fined as the ratio can represent the power at the true 
period relative to the overall level of noise present in the 
signal. The input signal is the same as for Fig. El but with 
the additional application of noise equal to half the am- 
plitude of the signal . In this case, the geometric sampling 
base used was x = 1.03, to minimise the effects of spectral 
leakage. The inset gives the amplitude of the correspond- 
ing window function, as a function of frequency; note the 
relatively minor interference contributions. 



to-noise normally quoted by astronomers, where the time- 
varying signal is superimposed upon a constant flux. For 
comparison, in a typical observing scenario where the full 
amplitude of the signal modulation might make up around 
10 per cent of the incident signal flux, our simulated noise 
would correspond to an observational combined signal-to- 
noise of around 33 and 10. 

In Figure ITU1 the noise level is comparable to the spec- 
tral leakage, and we see the optimal base choice clearly 
picked out by the A metric. For comparison, Figure 1111 
gives the RMS of the window function as a function of 
geometric base for the same dataset. Figure IT21 illustrates 
the degradation of the periodogram under the effects of 
extreme noise. Under these conditions, the choice of geo- 
metric base is much less important, because the amplitude 
of the noise is typically much greater than the amplitude 
of the spectral leakage. 

We conclude that noise is not a contributing factor 
to the choice of optimal base for standard astronomical 
applications. 

4.3. Randomised geometric spacing 

Given the favourable properties of a standard geometric 
sampling, we can ask whether the order of the sampling 
matters. We can preserve our choice of spacing while mod- 
ifying the order in which observations take place. This is 
equivalent to shuffling the gaps between observations. One 
motivation for doing this is that it allows individual ob- 
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Fig. 10. A for an N = 100 ordered geometric lightcurve 
with noise equal to 30% of the signal amplitude, as a func- 
tion of geometric base and injected signal frequency. 
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Fig. 11. RMS as a function of geometric base, for N = 
100. The optimum value is x = 1.032. 

servations to be placed with much more flexibility than 
an ordered geometric spacing, greatly easing scheduling 
problems. 

Figure EH plots the variation of RMS with geomet- 
ric base for a 25 point reshuffled geometric spacing. For 
comparison, the ordered geometric spacing is overplotted 
(dashes). In general, a randomised geometric spacing has 
a slightly higher RMS than the equivalent ordered spacing 
over the range of optimal bases. Figure PHI shows the win- 
dow function for the base with the lowest RMS. It should 
be noted that each random realisation has its own unique 
window function, and that for small values of N the RMS 
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Fig. 14. A typical window function for an optimum ran- 
domised geometric sampling of 25 data points, where 
x = 1.128. The RMS is 0.025. 



Fig. 12. A for an N = 100 ordered geometric lightcurve 
with noise equal to 100% of the signal amplitude, as a 
function of geometric base and injected signal frequency. 
High noise dominates the periodogram, washing out struc- 
ture in the plane. 
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Fig. 13. RMS as a function of randomised geometric base, 
for N = 25 (solid line). For comparison, the ordered geo- 
metric spacing has been overplotted (dashes). 

values can vary markedly from run to run. However, the 
range of variation in RMS across the space of possible 
random realisations shrinks rapidly with increasing N. 

The different behaviour between the ordered and ran- 
domised geometric spacing is more clearly illustrated in 
Figure ED The optimal range of bases remains almost the 
same as for the ordered case. In general, the smoothness 
(as measured by the RMS value) of the window function 
of the best randomised geometric spacing is not quite as 



good as the equivalent ordered spacing. However, the ran- 
domised spacing degrades much more gracefully with in- 
creasing geometric base - for sub-optimal choices of the 
geometric base, a randomised spacing outperforms the or- 
dered spacing. This has important implications for ob- 
serving programmes in which the total number of obser- 
vations is a highly dynamic quantity which cannot be ac- 
curately predicted in advance. In such a context, choosing 
a randomised geometric spacing would allow the observ- 
ing programme more flexibility with the total number of 
observations and their spacings, while seeking to optimise 
the schedule based on the current estimate of the run size 
and duration. Such an algorithm could form the basis of 
an adaptive dataset planning agent, an autonomous soft- 
ware entity which encapsulates sampling knowledge and 
attempts to exploit it in a real world environment. 

5. Practical application 

Figure presents the optimal geometric base for a range 
of randomly respaced geometric sampling scenarios. By 
taking the period range and run- length of interest, con- 
verting them into a relative limiting frequency, and then 
estimating the total number of observations to be made, 
the ideal geometric base for different observation regimes 
can be found. Expressing the limiting frequency in units 
of i^N, and substituting Equation ^ gives 



27V 2v r 



N 



N 



(10) 



As an example, an astronomer with access to a robotic 
network is planning a 3 week observing run, searching for 
periods ranging from 5 hours to 10 days. This corresponds 
to a minimum relative period of 0.0099, and thus a max- 
imum limiting relative frequency v re i of 101. If the total 
number of observations is likely to be around 100, then 
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Fig. 15. Randomised geometric base as a function of 
RMS, for N = 100 (solid line). For comparison, the or- 
dered geometric spacing has been overplotted (dashes). 
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Fig. 16. Optimal geometric bases, as determined by min- 
imum RMS, for datasets spanning 10 to 500 datapoints. 
For a relative maximum limiting frequency of 2 (lower 
panel), optimal bases tend to be closer to linear sampling 
than an equivalent sampling extending over a much higher 
dynamic range (upper panel). Simulations out to a relative 
maximum limiting frequency of 12.5 show almost no devia- 
tion from the optimal bases calculated for the z^nax/^N = 5 
case. 



(v/isn) ~ 2. Applying the lower curve in Figure ITfil the 
ideal sampling base is found to lie at around 1.02. 

Although FigurelTBlhas. in principle, to be recalculated 
for each particular value of (v/v^), in practice we find no 
change in the set of optimum bases for relative maximum 
limiting frequencies greater than 5, and the optimal min- 
ima are sufficiently broad that crude interpolation should 
suffice for smaller values. If the relative limiting frequency 



is found to be below 1, then enough datapoints exist to 
sample evenly - no further optimisation is required. 

Thus the observer should apply the following proce- 
dure to build an observing schedule. The signal range 
of interest must be identified, and the likely size of the 
dataset and its duration estimated. This allows the opti- 
mal geometric base to be calculated. Having ascertained 
the optimal base, Equation may be used to calculate a 
series of observation times. The resulting gaps may then 
be reordered as desired as the run proceeds, for example to 
accommodate practical constraints such as weather con- 
ditions or telescope maintenance, or to facilitate a single 
night's intensive coverage during a period of dark time. 

6. Conclusions 

We have found the best solution to date for placing lim- 
ited numbers of observations in such a way as to minimise 
aliasing effects arising from undersampling. By applying a 
simple geometric sampling rule, we are able to significantly 
outperform previous sampling schemes. Our scheme has 
the advantage that it is easy to apply at the telescope for 
datasets of arbitrary size, has calculated empirical solu- 
tions for different observing scenarios, and exhibits signif- 
icant flexibility arising from the invariance of the optimal 
sampling choice under the effects of reordering. As Figure 
[S] shows, by careful sampling the improvements that can 
be made in data quality, or alternatively, the savings that 
may be accrued in telescope time are substantial. 

Acknowledgements. The authors would like to thank Charles 
Williams for many useful discussions. ESS is funded through a 
PPARC e-Science studentship. This work is part of the eSTAR 
project, which supports AA and which is jointly funded by the 
DTI, PPARC and EPSRC. We are grateful to the Starlink 
project for their flexibility in seconding AA to eSTAR. We are 
also grateful to the European Union VOTech project. 

References 

Allan, A., Naylor, T., Steele, I., et al. 2004, in Proceedings 

of the SPIE, Vol. 5496, 313-321 
Bode, M. F., Carter, D., & Steele, I. A. 1999, Bulletin of 

the American Astronomical Society, 31, 1542 
Chaplin, W. J., Elsworth, Y., Howe, R., et al. 1996, 

Sol. Phys., 168, 1 
Deeming, T. J. 1975, Ap&SS, 36, 137 
Hessman, F. V., Gelderman, R., Naylor, T., Pennypacker, 

O, & Steele, I. 2004, American Astronomical Society 

Meeting Abstracts, 205, 
Home, J. H. & Baliunas, S. L. 1986, ApJ, 302, 757 
Littlefair, S., Naylor, T., Burningham, B., & Jeffries, R. D. 

2005, MNRAS, 358, 341 
Lomb, N. R. 1976, Ap&SS, 39, 447 

Mould, J. R., Hughes, S. M. G., Stetson, P. B., & Gibson, 

B. K. e. a. 2000, ApJ, 528, 655 

Nather, R. E., Winget, D. E., Clemens, J. C, Hansen, 

C. J., & Hine, B. P. 1990, ApJ, 361, 309 
Scargle, J. D. 1982, ApJ, 263, 835 



Eric S. Saunders et al.: Optimal placement of a limited number of observations for period searches 



9 



Schwarzenberg-Czerny, A. 2003, in ASP Conf. Ser. 292: 
Interplay of Periodic, Cyclic and Stochastic Variability 
in Selected Areas of the H-R Diagram, 383 — h 



